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Abstract. The Hardy-Littlewood majorant problem has a positive answer only for expo- 
nents p which are even integers, while there are counterexamples for all p ^ 2N. Montgomery 
conjectured that there exist counterexamples even among idempotent polynomials. This was 
proved recently by Mockenhaupt and Schlag with some four-term idempotents. 
(SJ " However, Mockenhaupt conjectured that even the classical 1 -I- e^'^'^ ± e^'^*('^+^)^ three- 

term character sums, should work for all 2fc < p < 2fc + 2 and for all /c G N. In two previous 
I papers we proved this conjecture for k = 0,1, 2, 3, 4, i.e. in the range < p < 10, p ^ 2N. 

Here we demonstrate that even the fc = 5 case holds true. 

Refinements in the technical features of our approach include use of total variation and 
integral mean estimates in error bounds for a certain fourth order quadrature. Our estimates 
I make good use of the special forms of functions we encounter: linear combinations of powers 

■ and powers of logarithms of absolute value squares of trigonometric polynomials of given 
degree. Thus the quadrature error estimates are less general, but we can find better constants 

p ^ '■ which are of practical use for us. 
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\0; 1- Introduction 

^ ■ Let T := R/Z. The Hardy-Littlewood majorization problem [6] is the question if for any 

CN ■ pair of functions /, (7 : T — t- C with |^| < / - that is, with / majorizing g - do we necessarily 

have II^IIp < 

Hardy and Littlewood noted that the Parseval identity easily implies this for all p G 2N 
^ ■ an even integer, but they also found that for p = 3 the property fails. Indeed, they took 

^ . / = 1 -|- ei -|- 63 and g = 1 ~ ei + (where efc(x) := e{kx) and e(t) := e^'^**) and calculated 

that II/II3 < WgWs- Later counterexamples were found by Boas [3] for all p ^ 2k and Bachelis 
[2J showed that not even allowing a constant factor Cp (i.e. requiring only ||5f||p < Cp||/||p) 
could save the property. 

Montgomery conjectured that the majorant property for p ^ 2N fails also if we restrict 
to idempotent majorants, see [lH p. 144]. (An integrable function is idempotent if its 
convolution square is itself: that is, if its Fourier coefficients are either or 1.) This has 
been recently proved by Mockenhaupt and Schlag in [10]. Their example is a four-term 
idempotent / and a signed version of it for g. For more details and explanations of methods 
and results see [3 [H] and the references therein. 
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In this paper we will be concerned with the even sharper conjecture, suggested by Mock- 
enhoupt in his habilitation thesis jSj- 

Conjecture 1. Let 2k < p < 2k + 2, where A; G N arbitrary. Then the three-term idempotent 
polynomial Pk '■= 1 + ei + ek+2 has smaller p-norm than Q/t := 1 + Ci — ek+2- 

Mockenhoupt presented an incomplete argument for the k = 1 case already in [9J. His 
argument hinted that some numerical analysis may be used in the proof, but we could not 
complete the solution along those lines. Nevertheless, we have proved this conjecture for 

= 0, 1, 2 in |E] and later even to = 3, 4 in [S]. 

One motivation for us was the recent paper of Bonami and Revesz [1], who used suitable 
idempotent polynomials as the base of their construction, via Riesz kernels, of highly con- 
centrated ones in U'iT) for any p > 0. These key idempotents of Bonami and Revesz had 
special properties, related closely to the Hardy-Littlewood majorant problem. For details 
we refer to [1] . For the history and relevance of this closely related problem of idempotent 
polynomial concentration in U see |H [5] , the detailed introduction of [7] , the survey paper 
py, and the references therein. The Bonami- Revesz construction, after suitable modification, 
directly and analytically gave the result for = 0. 

For larger fc, however, in (TJ [S] we used function calculus and support our analysis by 
numerical integration and error estimates where necessary. Naturally, these methods are 
getting computationally more and more involved when k is getting larger. "Brute force" 
numerical calculations still lead to convincing tables and graphes, but the increase of the 
number of nodes in any quadrature formula endanger the prevalence of theoretical error 
bounds due to the additional computational error, however small for reasonably controlled 
step numbers, but possibly accumulating for very large step numbers. 

Striving for a worst-case error bound incorporating also the computational error, we thus 
settled with the goal of keeping any numerical integration, i.e quadrature, under the step 
number N = 500, that is step size h = 0.001. Calculation of trigonometrical and exponential 
functions, as well as powers and logarithms, when within the numerical stability range of 
these functions (that is, when the variables of taking negative powers or logarithms is well 
separated from zero) are done by mathematical function subroutines of usual Microsoft 
Excel spreadsheet, which computes the mathematical functions with 15 significant digits of 
precision. Although we do not detail the estimates of the computational error of applying 
spreadsheets and functions from Microsoft Excel tables, it is clear that under this step 
number size our calculations are reliable well within the error bounds. For a more detailed 
error analysis of that sort, which similarly applies here, too, see our previous work [7j, in 
particular footnote 3 on page 141 and the discussion around formula (22), and see also the 
comments in the introduction of [8]. 

We keep using the fourth order quadrature formula, presented and explained in |8], see [H 
Lemma 5]. However, another new argument also has to be invoked for k = 5 compared to 
k = 3,4, because in this case the anal5d;ic scheme of proving fixed signs of certain derivatives 
simply break down. Using the special form of our integrands and the resulting form of 
estimates with the initial trigonometrical functions, we thus invoke the special quadrature 
error estimate of Lemma [71 These estimates make good use of the concrete form, local 
maximum values and alike, of the functions G* in question, but the theoretical estimates 
with Var G (the total variation of the function G) and ^ G^iO over local maximum values 
( oi G, might have some theoretical interest, too. 

Second, as already suggested in the conclusion of [7] and applied in f8] for = 4, we 
use Taylor series expansion at more points than just at the midpoint to := k + 1/2 of the 
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t-interval {k,k + l), thus reducing the size of powers of (t — to)) from powers of 1/2 to powers 
of smaller radii. 

Finally, we needed a further consideration in proving that the approximate Taylor poly- 
nomial P{t), minus the allowed worst case error 6, still stays positive in the interval of our 
Taylor expansion. Basically, in we could always use that the polynomials p{t) := P{t) — 6 
were totally monotone - now some occurring approximate Taylor polynomials will not have 
this feature, and we need a more refined calculus to succeed in proving their constant sign 
over the interval of investigation. 

Key to this is the consideration of the variance of some of the derivatives of p, for if a 
function vanishes somewhere inside an interval, than its variance exceeds the sum of the ab- 
solute values taken at the left and right endpoints of the interval considered. This elementary 
fact comes to our help in concluding that p{t), and hence the considered difference function 
d{t), approximated by P{t) within a certain error 6, keeps constant size; for if the first j 
derivatives are positive at the left endpoint, and the jth derivative preserves the positive 
sign all over the whole interval, then there is no way for p to vanish anywhere in the interval. 

2. Boundary cases of Conjecture [T] at p = 2fc and p = 2k + 2 
Let A; G N be fixed. (Actually, later we will work with k = 5 only.) 

We now write F±{x) := l + e(x) ±e((A; + 2)a;) and consider the p^^ power integrals f±{p) := 
Jq \F±{x)\^dx as well as their difference A(p) := f_{p) — fj^{p) := \F_ix)\'^ — \F^{x)\Pdx. 
Our goal is to prove Conjecture [H that is A{p) > for all p G {2k, 2k + 2). 

Let us introduce a few further notations. We will write t := p/2 G [k,k + 1] and put 

1/2 

G±(x) := |F±(x)p, g^{t) := ^/±(2t) = j Gi{x)dx, 



1/2 

(2) d{t) := ^A{2t) = g.{t) - g+{t) = J [G'_ix) - G^x)] dx. 



Formula also yields that denoting Htj ±{x) := G*j_{x)log' G±{x) the explicit integral 
formula 

S^\t) = g^l\t) - g^i\t) 



(3) 

holds true, and so in particular 

Li [0,1/2] + II -f^tj-IUi [0,1/2] (j G N). 

We are to prove that d(t) > ior k < t < k + 1. First we show at the endpoints d vanishes; 
and, for later use, we also compute some higher order integrals of G±. Actually, here we can 
make use of the following lemma, already proven in [HI Lemma 3]. 



1 



1/2 1/2 

Giix) log-'' G-{x)dx - [ G*+(x) W G+ix)dx 



1/2 1/2 

J Htj-{x)dx — J Htj^+{x)dx. 
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Lemma 2. Let p E N with 1 < p < k + 1. Then we have 

p-(k+2) 

^ a±(z/)e^ 



(5) 
where p : 

(6) 



G<i 



u=0 



with 



±11 



p- p 
A 



k + 2 



and A := u — p{k + 2) is the reduced residue of u mod k + 2. Therefore, 



1/2 



p-(k+2) 



\G±\ 



u=0 



In particular, J^^^ \G+\^ = j^"^ for all < p < k + 1 and thus d{k) = d{k + 1) = 0. 

Apart from the immediate result that d vanishes at the endpoints of the critical interval 
[fc, /c + 1], we will make further use of the above explicit computation of p^^ power integrals 
of G. To that we need the precise values of these square sums of coefficients, which is easy 
to bring into a more suitable form for direct calculation. Namely we have 

p-{k+2) 



'1/2 



A{p) :- 



p- p 

A 



/i=0 



2 P-P / \ 2 

p- p 



A=0 



A 



At=0 



2p-2p 



[k+2\ 



A~i/-M(fc+2) 

(7) = 1,3,15,93,639,4653,35169 for p = 0, 1, 2, 3, 4, 5, 6, respectively. 

Corollary 3. For all p < (/c + 1) we have jj^^^ G^. = \A{p) with the constants A[p) in ([7]). 

With the aid of these explicit values, even arbitrary power integrals of G± can be estimated. 
Proposition 4. Let p G N and p <k + 1. Then with the constants A{p) in ([7]) we have 

1/2 1/2 

(8) j Gl<\9^-^A{p) {r<p) and j GI<\a^Ip{p) {r > p). 



Proof. As < G < 9, for the first estimate one can use < 9^^^C. The second estimate 
is directly furnished by Holder's inequality with exponents p = p/r > 1 and g = l — 1/p. □ 

3. Analysis of G± 

To start the analysis of G{x) := G^^±{x), let us compute its ^-derivatives. As in formula 
(7) and the following lines of [8], in case = 5 we find easily 

(9) G±{x) = 3 + 2{cos(27rx) ± cos(127rx) ± cos(147rx)} 

G^l'^^^\x) = (-4)'"+i7r2™+i . {sin(27ra;) ± 6^™+^ sin(127rx) ± 7^""+^ sin(147rx)} , 
G^i'^\x) = 2(-4)™7r2'" . {cos(27rx) ± 62"^ cos(127rx) ± 7^™ cos(147rx)} . 
Consequently we have 



(10) ||G±|U<9=:M( 



0; 



\G 



(m) I 



< 2"^+V"{l + 6*" + 7™} =: M„ 



1,2,...), 



that is IIG 



(m) I 



<Mm{m = 0,1,2,3,4) with 



;il) Mo = 9, Ml = 175.929... < 176, Ms = 6790.287... < 6800, 

M3 = 277, 816.239... < 280, 000, M4 = 11 527, 002.2... < 11, 600, 000. 
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Lemma 5. Both functions and G^{x) have seven local maxima in T = (— |, |]. At 

zero-symmetric pairs of maximum places of these even functions the same maximum values 
occur, so presenting these values with multiplicity 2, they are the following: 



G^{x) 


G^x) 


C 




multiplicity 


c 


G-iO 


multiplicity 





9 


1 


^ ±0.076 


< 8.662 


2 


^ ±0.151 


< 7.701 


2 


^ ±0.227 


< 6.279 


2 


^ ±0.302 


< 4.628 


2 


^ ±0.377 


< 3.005 


2 


^ ±0.448 


< 1.661 


2 


0.5 


1 


1 



Proof. As G is a degree 7 trigonometric polynomial, if it has n local maximums, then there 
are the same number of interlacing minimums, so altogether 2r;, < 2 degG' = 14 roots of G'. 
So the number of maxima is at most 7, which we will find - taking into account evenness of 
G, and thus the same symmetrically located maxima and minima in [—1/2, 0] and in [0, 1/2] 
- so no further local maxima can exists. 

As G is even, it suffices to analyze [0, 1/2]. We start examining the functions by tabulating 
it with step size h = 0.001, and identifying the indices i where monotonicity of the G{xi) 
turns from increase to decrease. Then there has to be a local maximum at some point 
Ci G [xi-i,Xi+i]. Clearly one of the nodes x- G Xj, Xj+i} has |Cj — < h/2. The second 

order Taylor expansion around Q now gives G{(i) — G{x'^ < || |oo(|)^, as G'{Ci) = 0. 
Applying (IIO]) Ma < 6800 and h = 0.001, we obtain G(Ci) - G{x'i) < 0.00085 < 6 := 0.001. 

In the table above we recorded d ~ Xi with error < h = 10~^ and an upper estimation of 
the corresponding maxima using < G'(x^) + 5 < max G{xi), +S. □ 

Denote Var (■?/', [a,b]) the total variation of the function ip on [a,b], and in particular let 
Vai{tp) := Var('?/',T). As an immediate corollary to the above lemma, we formulate here 

Corollary 6. Denote by Z := Z± the set of local maximum points of G = G±. For any 

positive parameter t > Q we have Var(G*) < 2 ^^^^^ In particular, Var(G±) < 74. 

Proof. It is easy to see that for a piecewise monotonic function ip one has Var('?/', [a, 6]) = 
Var(|V|, [a, 6]). It follows that Var(|V^|, [a, 6]) = |^'|. 

Furthermore, for a piecewise monotonic function, like G or G*, the total variation is the 
sum of the change of the function on each of its monotonicity intervals. Since we are talking 
about periodic functions, i.e. functions on T, with only finitely many critical points, it is 
clear that the local maximum and minimum places - with the latter denoted by f2 C T, 
say - interlace and monotonicity segments connect these neighboring local extremum places. 
Therefore the total sum of all the changes is 

Var(G*) = - = 2 E " 2 E ^*(^) < 2 E 

e,rjezun c&z weO C&Z 

(e,r7)n(zuQ)=0 

taking into account G > 0, too. 

Whence the first assertion of the Corollary, while the last is just a small calculation adding 
the maxima (taken into account according to multiplicity) in the columns of the table of 
maxima in Lemma O □ 

4. Estimates of \H{x)\ and of HiJ^^Hoo 



Let us start analyzing the functions 
(12) H{x) := ift,,,±(x) := G\x) \og' G{x) (x G [0, 1/2]) (t g [A;, A; ± 1], j G N). 
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To find the maximum norm of Htj^j-, we in fact look for the maximum of an expression of 
the form f*| logf p, where v = G{x) ranges from zero (or, if G 7^ 0, from some positive lower 
bound) up to ||G||oo < 9. 

A direct calculus provides a description of the behavior of the function a{v) := as,m{v) : = 
v^l logf p for any s > and m G N on any finite interval [a,b] C [0, 00), see 0, Lemma 6]. 

For the application of the above quadrature ffTSjl we calculated (c.f. also [HI, (15)]) 

H"{x) := H'lj^^{x) = G"{x)G'-\x) log^"^G(x) {tlogG(x) + j} 

(13) + G'\x)G'-'{x) \og^-' Gix) {tit - 1) log^ Gix) + j{2t - 1) log G{x) + j(j - 1) } • 

However, the error estimation in the above explained quadrature approach forces us to 
consider even fourth x-derivatives of H = Htj^± using H^^ = ^^=o (m) G)'-'^"™-'. 
We have already computed in fS] respective formulae for (G*)^™'^ and (logG-')^™^ for m = 
1,2,3,4 (c.f. [SI (17), (18)]). Substituting these in H^^ resulted in the general formula [HI 
(19)] stating with L := logG 

H'v ^ g'-'G''{jU - l)(j - 2)(j - 3)U-' 

+ [4t - 6]j(j - l)(j - 2)L^-' + [6t^ - 18t + ll]j(j - 
+ [2t^ - 9t^ + lit - 3]2jL^-^ + t{t - l)(t - 2){t - 3)U^ 

+ 6 ■ G'-'G"G"{jU - l)(j - 2)L^--' + 3{t - - 1)L^-' 

(14) + [3t^ -6t + 2)]jL^-^ + t{t - l)(t - 2)L^} + G'-^G^^ {jL^'^ + tL'} 
(3 ■ G'-^G"^ + 4 ■ G*-2G'G'") |j(j - l)L'-^ + {2t - + t{t - 1)L^}. 

Finally, from that and writing in ||G*^'^^||oo < we were led to the general estimate 
[F, (20)]. As now the values of Mm are estimated by (fTTj) . the corresponding values can be 
written in, and putting also £ := \L\ = \ logG| formula [8, (20)] yields 

IH'^^I < 959, 512, 576 ■ v'-'[j{j - l)(j - 2)(j - 3)&-' 

+ [4t - 6]j(j - l)(j - 2)F-3 + [6^2 _ I8t + ll]j(j - 
+ [2t^ - 9t^ + lit - 3]2jf-^ + tit -l){t- 2){t - 3)f| 

(15) + 1, 263, 820, 800 ■ v'-^'^jij - l)(j - 2)f-^ + 3{t - l)j{j - l)f~' 

+ [3t2 - 6t + 2]jf-^ + t{t - l)(t - 2)P^ + 11, 600, 000 • v'-^ {jP-^ + t&} 
+ 335, 840, 000 ■ v'-^\^j{j - 1)£-''"^ + (2t - 1)3 f'^ + t{t - 1)^}. 

On the other hand, for reasons becoming apparent only later from the improved quadrature 
error estimate in Section [5l here we need to derive another consequence of formula f|T^ . We 
now substitute the norm estimates of ffTTl) by M^'s into f|T^ only partially, that is, we leave 
(apart from all powers of G) even one of G' without estimation by Mi, wherever G' occurs, 
in order to take advantage of our quadrature utilizing expressions of the form G^IG'log"* G|. 



Inserting k 
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5 and the numerical values of Mi , M2 , M3 and M4 from fill I) this leads to 



\H 



IV \ 



< 5,451, 776 ■ G'-^\G'\[jij - l)(j - 2)(j - 3)f~^ 

+ [At - Q]j{j - l)(j - 2)F-3 + [Qt^ - 18t + ll]j(j 



1)F~ 



+ [2r -9r + llt- 3]2j&-' + t{t -l){t- 2){t - 3)& 



+ 7, 180,800 ■ G'--'\G'\\j{j - l)(j - 2)f-' + 3(t - l)j(j - 



+ [3t' - 6t + 2)]jf'^ + t{t -l){t- 2)f 



+ 1, 120, 000 ■ G'-^\G'\\j{j - + {2t - l)jf-^ + t{t - l)& 



logGI). 



(16) 



+ 138, 720, 000 ■ - l)^'^ + (2t - + t{t - l)& 

+ 11, 600, 000 ■ G^-^ + tP] (with, as always, 

5. Quadrature with variation 

In the paper [7] we used Riemann sums when numerically integrating the functions H := 
GMog-'G along the x values. A new feature of the subsequent paper [8], among other 
things, was the application of a higher order quadrature formula- Namely, in [8], formula 
(12) and (13) we recalled the following easy-to prove elementary fact. Let be a four times 

continuously differentiable function on [0, 1/2], G N, h := 1/(2A^) and denote x„ : = 
for n = 1, 2, . . . , A^. Then we have 
1/2 



2n 



AN 



(17) 



n n=l ^ 



ip" (x^ 



TV 



< 



60 ■ 210 A^5 



max 

\^-^n\<2 



^^^(a;)|. 



In [8] we then used this with the further obvious estimate maX|^_^,^|<h |y) (x)| < \\if 

II /Vll 

resulting in the further estimation of fflTI) by go.210 • 
1/2 



i n=l 



2n 



AN / 2A^ 



2n - 1 



AN / 192A^3 



< 



60 ■ 210 A^4- 



We intend to use the quadrature formula (fT8|) to compute approximate values of (i'(5), c?"(5), 
d"'{5) and then even rf(^)(to) with various values of to ^ [5, 6] and j G N. However, use of 
direct estimations of ||-f^/(,j-,±||oo in the quadrature would result in step numbers as high as 
800, already inconveniently large for our purposes. Thus here we invoke a further, more 
detailed analysis of the quadrature formula, aiming at bounding the step number further 
down below 500 with the improved error estimation. 

The basic idea is that we try to apply (1171) directly. For continuous f^^, the local maximum 
are attained at certain points ^„ G [x„ — h/2,Xn + h/2], and the error bound becomes 
Sn=i lv^^^(^")l- f^^^ ^his sum is a Riemann approximate sum of the integral (and not the 
maximum) of the function ip^'^ , so we will get approximately 2A^ ■ J 



V2 ,^IV 



That is, we 



arrive at the norm, instead of the L°° norm, of the function 'p^'^ . 

So we try to make use of this observation for HtQj^± in place of ip. Again, direct estimation 

of the error in this approximation J2n=i l-^/oJ,±(^"^)l ~ '^■^ ' lo^"^ l-^/oj',±l' ^^^^ theoretically 
possible, does not provide nice and numerically advantageous results. Instead, we estimate 
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the function H^^j_^, similarly as above, with functions involving G, logG and even deriva- 
tives of G, and then spht the estimation of the sum "^^^i |i?/(,j _|_(^n)| to estimations of 
similar Riemann sums of such simpler functions. For such combinations as log"* G{x) 

or G^{x) log-' G{x)G'{x), we will find suitable error bounds and explicit computations or esti- 
mations of the L^-norms, finally resulting improved estimations of the error in the quadrature 
formula. More precisely, we can derive the following improved special quadrature estimation, 
which subsequently will be used for -ffto,j,± (with various j and to) in place of ip. 

Lemma 7. Let Br, Dr > 0,1 < tr < T and jr > for r = 0,1, . . . , R. Assume 

R 

(19) |/^(x)| < J2 [BrC'-ix) ■ I logG(x)p'- + DrG'^{x)\G'{x)\ ■ I logG'(x)|^'-} . 



r=0 



Then for arbitrary N eN the quadrature formula 



1/2 



LP 



(20) 

holds true with 



N 

E 

n=l 



2n- 1 



4N / 2N 



2n - 1 



1 



AN / 192A^3 



< 



1 



R 



60 • ATS 



{BrQN{G,tr,jr) + Dr Q*n{G ,tr , jr)} 



r=l 



(21) Qn{G, t,j) := xU ^ 0) (^max v'\ \ogv\^^ N + \og^ 9 | iV ^ G* + ^ Var(G* 



and 



Q*^{G,t,j) := xU ^ 0) ( max t;*| logw| 

y [u,i/ yj 



14 



+ 1700 



(22) 



+ log^' 9 



N 

TTT 



Var(G*+^) + 88 Var(G*) + 1700 



Q2t 



\ T 



Proof. As in the preceding arguments, we denote x„ := {2n — 1)/(4A^) and h := 1/(2A^) for 
n = 1, . . . , and even for n = 1 — A^, . . . , A^. By the inequality (I17p . the condition (IT^ and 
making use that all the arising terms in this estimate are continuous and even, we find with 
some appropriate, symmetrically chosen ^„ G [x„ — h/2,Xn + h/2] that 

1/2 

h 

< — — — max \ip^'^{x) 



I ^ ~^\v{xn)h + ^"{xn)^^ 



60 ■ 25 



< 



N 



^ , max 
60 ■ 25 ^ I i<h 

n=i I"" ^"1-2 



R 



J2 {BrG'^' (x) I log G{x) + DrG'^ {x)\G'{x) 1 1 log G{x) } 



r=0 



(23) 



1 



N R 



60 ■ 25 2 



J2 E i^rG'^ (^-) I ^(^-) I'" + ^rG'^ i^n)\G'{U) \ \ log G{^n) |> } 



AT- 
60^2 



n=l-N r=Q 
R 



N 



N 



uJ2\Br G*^(en)|l0gG(Ur + E G'^{U\G'{U\\l0gG{^r, 



\3r 



r=0 



n=l-Af 



n=l-Af 
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So we are left with the estimation of the inner sums. There are two type of sums here, 
the first being without and the second with its appearance. For a more concise 

notation let us introduce the exponent k e {0, 1}, and then consider the generic inner sum 

S Sit,j,K) := En=l-7vG'*(en)|G'(en)|1l0gG(en)K. 

To start with, when j = we can directly compare this sum to the corresponding integral. 
Recall that for any function ip of bounded total variation Var(-?/') := Var(-(/', [a,b]) on an 
interval [a, b], and for any partition of [a,b] a.s a = Xq < Xi < ■ ■ ■ < Xi < ■ ■ ■ < xm-i < xm = b 
with the fineness of the partition 5 :— maxi=i^,,,^M{xi — Xi-i) and with any selection of nodes 

9i e [a;i_i,Xi], the Riemann sum ^fLii'{^i){xi ~ Xi-i) approximates '^{x)dx within the 
error 5Var('0). So we obtain 

N 

(24) S{t,Q,K)= ^ G\^n)\G'{^nW <2N G^\G'\^ + 2NhYw{G*\G'\^). 

n=l-N 



If K = 0, then the first term is 2N G*, and for k = 1 it is nothing else than 2N Vai^j^G^^^) 
on T. As 2Nh = 1, for /t = the second term is Var(G*), while for k = 1 we can also obtain a 
similar type estimate using that Var(|\E'|, [a, b]) = Var(^', [a, b]) = |^''|. Namely we obtain 

Var(G*|G'|) = J |(G*G')'| < j tG'-^G'^ + j G'\G"\ < Mi j tG'-^\G'\ + j G'\G"\ 



j G2* j G"2 < 176 Var(G*) + 3400 ^ 



(25) <176Var(G*) + 

with an application of the Cauchy-Schwartz inequality and computing 



Q2t 




9 /I + 362 + 492 
Stt^ • \ ^ 3395.144... < 3400. 



So collecting terms furnishes 

J 2N G* + Var(G*) if k = 0, 

(26) S{t, 0, k) < |27V_1_ Var(G*+i) + 176 Var(G*) + 3400^/^ G2* if At = 1. 

Observe that the right hand side of this estimate is just 2Qn{G, t, 0) and 2Q*j^{G, t, 0) when 
At = and 1, respectively, so the part of the assertion for j = is proved. 

For J > we estimate S{t,j, k) by first cutting the sum into parts according to ^„ G X : = 
{x e T : < G{x) < 1/9} and ^„ ^ X. The first of these sums can then be estimated 
by maxo<^<i/9 v*| log'' v\ ■ X^^^gx l^'(Cn)r, the sum being < constant 2N for At = while for 
At = 1 approximately 2A'" /^^ |G'| = 2A'" j^\G'\xx, where xx is the characteristic function of 
X. 

That latter integral of |G'| on X is just the total variation of G{t) along its segments 
of range between and 1/9. More precisely, as G is a trigonometric polynomial, hence 
piecewise smooth with at most (actually, exactly) 2degG = 14 monotonicity intervals Im 
(m = 1, . . . , 14) within T = U^^^/m, this whole total variation can amount at most 14 times 
the maximal possible variation from to 1/9 on each part of X belonging to one monotonic 
segment That is, |G'| = Ei!=i /xn/^ 1^1 = EL=i Var(G, X n J^) < 14 • 1/9. In all, 
the contribution of the main term 2N \G'\ is at most 14/9 • 2A'". (In reality, that variation 
is numerically even less, but this term will not be too interesting anyway.) 
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Next we apply the general Riemann sum error estimate to to infer J2^„ex — 

En=i-N\G'{in)\xx{U<2Nj^\G'\xx + 2N-h-VaT{\G'\xx) < yiV + Var(G'xx). 

We now show that this latter variance does not exceed Var(G'). In view of the additivity 
of the total variation on intervals, Var(G'xx) = J2m=i'^^^iG'xx , Im) , so it suffices to prove 
Va.T{G'xx, Im) < Var(G",/m). Recall that Im = [am,bm] is, by construction, one of the 
intervals of monotonicity of G, hence a segment of T where G' has constant sign, with zeroes 
(and sign changes) of G' at both endpoints. Then either G{am) is a local minimum of G 
and G{bm) is a local maximum of it, or conversely, corresponding to the cases when on Im 
G' > or G' < 0, respectively. By symmetry, we can restrict to the first case, when G is 
increasing on 7^. If G > 1/9 on /„, that is, if already G{am) > 1/9, then n X = 
and Vai^G'xx, Im) = < Var(G",/m). Also if G < 1/9 on the whole interval /„, then 
Vai^G'xx, Im) = Var(G', /„). The only case when n X is nontrivial is when fl X = 
[am,Cm] with ttm < Cm < &m and G{am) < G{cm) = 1/9 < Gipm)- In this case, however, 
G'xx = G' on [am,Cm[, has a jump from G\cm) to at c^, and constant zero afterwards 
until the end of the interval 7^, whence Var(G'xx, Im) = Var(G', [a^, c^]) + \G\cm) — 0| = 
Var(G', K,c^]) + \G\cm) - G'{hm)\ < Var(G', [a„,c„]) + Var(G', [c^,M) = Var(G",/™). 



So indeed we have Var(|G"|xx) < Var(G'). Finally, using the above estimation of J J G""^, 



Var(IG'lxx) < Var(G') = J \G"\ G"^ < 3400. 

Writing in the maximum of v*'\ logf p, collection of terms results in 

(27) y G\in)\G'{in)U logG(eJP' < i max v'\ \ogvA \ 1^ if = 0, 
^ ' f-'^ ' ^ - VI0.V9] ' ';lfAr + 3400 ifK = l. 

In the second sum over ^„ ^ X we have G(^„) G [1/9,9], hence |log(j>'(^„)| < log 9, so 
bringing out this estimate from the sum and then extending the summation to all n leads to 

N 

(28) 5^G*(a)|G'(a)r|logG(a)K<log^'9 ^ G*(a)|G'(a)r = W 9 5(t, 0, «:). 

i^iX n=l-N 

Summing up, if j 7^ then the upper estimate of 

N 

E G\U\G'{Un^ogG{U\' 



n=l-N 



(29) 




max[o,i/9] log wp) 2X + log^' 9 ■ 2Qn{G, t, 0) if k = 0, 

(max[o,i/9]t^*|logt;p) {f X + 3400} + log^' 9 • 2g^(G, t, 0) if k = 1. 



follows for the generic term, and so taking into account the notations (l2Ti) and (l22l) . from 
(1291) and (|23D the lemma follows. □ 



6. Derivatives of the difference function d{t) at the left endpoint 
With the improved quadrature we now calculate the values of d'{5), d"{5) and d'"{5) first. 
Lemma 8. We have d'(5) > 0. 

Remark 9. Actually, d'{5) = 0.00287849... by numerical calculation, but formally we don't 
need an a priori knowledge of the value. Of course, putting together the argument we needed 
to take it into account, but the proof of the Lemma is deductive. 
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Proof. From ( TTB]) . substituting j = 1 and t = 5 and denoting, as elsewhere i := \L\ = \ \ogG\ 



(30) 



W^\< G|G"|{839, 573, 504 + 654, 213, 120£} + G^|G"|{337, 497, 600 + 430, 848, 000£} 
+ G^\G'\{10, 080, 000 + 22, 400, 000£} + ^^{1, 248, 480, 000 + 2, 774, 400, 000£} 
+ {11, 600, 000 + 58, 000, 000£} . 



This estimate is of the form of condition f ll9p . suitable for the application of our improved 
quadrature in Lemma [TJ which we invoke with := 500 here. Therefore, we compute the 
expressions f l2T]) and f l22|) with = 500 and with the occurring pairs of values of t and 
j = 0, 1 as follows. 

First of all, observe that according to [8\ Lemma 6] for j = 1 and t = 1,2,3,4 we have 
max v*| logf p = 9~*log-' 9 (as the maximum place vq = exp(— j'/t) is larger, than 1/9). For 

the computation of G* and Var(G*) we refer to Corollaries |3] and El These lead to 
Qlooi^, 2, 0) 



Qlooi^^ 2, 1 
Qlooi^y 3, 
Qlooi^' 3, 1 
QbooiG, 3, 

Q50o{G, 3, 1 
^500(^,4,0 

Q50o(G,4,l 



< 123, 987 + 6509 + 6585 = 137, 081, 



< ( + 1700 1 + log 9 ■ 137, 081 < 301, 803, 
9 \ 9 



< 616, 734.71 + 43, 643.12 + 42, 973.37 = 703, 352, 



< ( + 1700 I + log 9 ■ 703, 352 < 1, 545, 490, 
9"^ \ 9 / 

< 3, 632, 988 + 325, 636 + 318, 808 = 4, 277, 432, 

< ( + 1700 I + log 9 ■ 4, 277, 432 < 9, 398, 487, 
9'^ \ 9 / 

< 46,500 + 1851 = 48,351, 
log 9^ 



93 



-500 + log 9 ■ 48, 351 < 106, 240, 



< 319, 500 + 14, 532 = 334, 032, 
log 9^ 



< 



94 



-500 + log9 ■ 334, 032 < 733, 944. 



It remains to apply Lemma [7] both for iJ+ and H_ with the coefficients Br,Dr read from 
f l5U]) and the corresponding Q5oo{G,t, j), QlQQ{G,t, j) estimated according to the above list. 
Executing the numerical computations leads to 



(31) 



1/2 



500 

E 

n=l 



H. 



± 



2n 



1 



2000 / 1000 



2n - 1 



1 



2000 / 192 ■ 5003 



< 0.0009745.. 



Thus the quadrature approximation to integrals of H± lead to approximate values within 
the error 6 := 0.001. This error estimation is applied to both and H . The approximate 



value ofrf'(5) = /^ ^ iJ_ - /^ ^ i7+ from the quadrature is found to be 0.002878492... > 0.002, 
while the total error incurred is still bounded by 26. Therefore, d'{5) > 0.002 — 26 = and 
the assertion is proved. □ 

Lemma 10. We have d"{5) > 0. 

Remark 11. By numerical calculation, d"{5) ~ 0.033815603. 
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Proof. Now we want to use the improved quadrature again, hence we start with substituting 
t = 5, j = 2 into formula (ITB]) to derive 



(32) 



< G|G"|{774, 152, 192 + 1, 679, 147, 008£ + 654, 213, 120r} 
+ G^\G'\{172, 339, 200 + 674, 995, 200£ + 430, 848, OOO^^} 
+ G^|G"|{2, 240, 000 + 20, 160, 000£ + 22, 400, OOO^^} 
+ G^{277, 440, 000 + 2, 496, 960, 000£ + 2, 774, 400, OOO^^} 
+ {23, 200, 000£ + 58, 000, OOOf} . 



Now we may set the step number to = 400. The values of the occurring Q^ooiG, t,j) and 
Qloo{G,t, j) can now be estimated as follows. 



Q:oo(G',1,0)<112,282 
Q:oo(G',2,0) < 580,005 
QlooiG, 3,0) < 3,550,835 
g4oo(G',3,0) < 39, 051 
Qioo{G, 4, 0) does not occur 



g:oo(G',l,l)< 247,274 
g:oo(G,2,l)< 1,274,463 
g:oo(<^,3,l) < 7, 801,987 
^400(^,3,1) < 85, 804 
Q4oo(G',4,l) < 593,541 



Q:oo(G,1,2) < 543, 316 
Qloo{G,2,2) < 2,800,281 
Q:oo(G',3,2)<17, 142,718 

g4oo(G',3,2) < 188,530 
Q4oo(G',4,2) < 1,304,143 



Applying Lemma [7] for either or iJ_ with the coefficients -B^, read from 
corresponding Q^{G,t, j), Q*j^{G,t,j) above, the numerical computations yield 



and the 



(33) 



1/2 



400 



n=l 



2n - 1 



1 



4 ■ 400 / 2 ■ 400 



HI 



2n - 1 



1 



4 ■ 400 / 192 ■ 4003 



< 0.0071... =: 5. 



This quadrature error estimation is applied for both if+ and so the total error incurred 
is still bounded by 25, while the approximate value of d"{b) = J^^^ if_ — f^^'^ from the 
quadrature is found to be 0.033815603. Therefore, d"{5) > 0.033815603 - 25 > 0. □ 

Lemma 12. We have d"'{5) > 0. 

Remark 13. By numerieal caleulation, d'"{5) ~ 0.183547634.... 

Proof. Now it suffices to apply the less refined estimates from f|T5|) with t = 5, j = 3 to get 

\H^^{x)\ < 959, 512, 576 • f {84 + 426£ + 462^^ + 120i^} 
(34) + 1, 263, 820, 800 ■ v^{6 + 72i + Ulf + 60f} + 11, 600, 000 ■ {3f + 5i^} 

+ 335, 840, 000 ■ v^{6£ + 27i^ + 20f}. 

In this estimation all the occurring functions of type v'^i"^ have maximum on [0, 9] at the 
right endpoint v = 9 in view of [SI Lemma 6]. Therefore we can further estimate substituting 
£ = log9 and v = 9. Thus we finally obtain \H^^{x)\ < 2.82932 • 10^^. 

To bring the error below 6 = 0.091 we chose the step number N large enough to have 



14 



2.83-10 
60 • 210 



< 6 



I.e. 



iV > A^n 



2.83 ■ 1014 
60 ■ 210 ■ 0.091 



475... 



Calculating the quadrature formula with = 500, we obtain the approximate value d"'{5) 



/o - /q'^^ H+ ^ 0.18354763424..., whence d"'{5) > 0.18354763424... - 2 • 0.091 > 0. □ 



'1/2 
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7. Signs of derivatives of d{t) and conclusion of the proof of Conjecture [T] 



After examining the values of derivatives of d at the left endpoint t = 5, now we divide 
the interval [5,6] to 3 parts. First we will prove in Lemma ITSl that d^^{t) > in [5,5.13]. 
In view of the above proven Lemmas El [10] and [121 it follows, that in this interval also 
d"'{t),d"{t),d'{t) > 0. 

Next we will consider d' in the interval [5.13, 5.72]. Lemmas [T7I and [T9] will furnish d' > 
also in this domain. Consequently, d is increasing all along [5, 5.72], and as d{5) = 0, it will 
be positive in (5,5.72]. Finally we will show Lemma f2U\ giving that d{t) is concave in the 
interval [5.72, 6]. As d{5.72) > and d{6) = 0, this entails that the function remains positive 
on [5.72, 6), too, whence d > on the whole of (5, 6). 

Now we compute a good approximation of d'^{t) on the interval [5, 5.13] and using it show 
that d'^{t) stays positive in this interval. 

In [5,5.13] the fourth derivative of d{t) has the Taylor approximation 

(35) d'^{t) = ~ + ^"(^^^' 5-065, t), where 

j=o ^■ 

Rn{d'^, 5.065, t) := ^ (t - 5.065)"+' . 
Therefore using we can write 

|i?„(rf^^,5.065,t)| < + l|g^.n+5,-|UM0.1/2l 

I nv > ;i - (n + 1)! 

(36) < kn,n+,M\oo + kn,n+5,-\\oo q ^q^u-,! 

^ ^ ~ {n + l)\ 

^ maXj^_5,065|<0.065 ||-f^5,n+5,+ ||oo + maX|g_5.065|<0.065 || -f^f,n+5,- 1| oo ^ 065"+' 

~ (n + iy. 

So once again we need to maximize ([T2!) . that is functions of the type v^ \ logf I*", on [0,9]. 
From IHl Lemma 6] we get, say for all n < 30 

(37) max WH. ^+5 ±(x)\\^ < max max ^;«| logt;r+^ = 9^+^ log"+^ 9. 

5<C<5.13 ^' • ^ 5G[5,5.13] i)e[0,9] 

Choosing n = 6 yields ]|if^^„+5^±(x)]|oo < 452, 775, 589, and the Lagrange remainder term ([36]) 
of the Taylor formula §^ can be estimated as \Rn{d^^,t)\ < 0.0008808... < 0.0009 =: 67. 

Now we have to calculate the value of d^^^ (t) - that is, the two integrals in (j3|) - numerically 
for k = 5, t = 5 and j = 4, 5, . . . , 10 to determine the Taylor coefficients in the above 
expansion. However, this cannot be done precisely, due to the necessity of some numerical 
integration in the calculation of the two integrals in formula ([3|). We apply our numerical 
quadrature to derive at least a good approximation. 

Denote dj ^ dj^4 the numerical quadrature approximations. We set 6 := 0.187 and want 
that \d^^{t) - Peit)] < 5 for 

(38) Pg(t):=^^(^_5.065)^ 

j=0 ^' 

In order to achieve this, we set the partial errors 5q,. . .,5q with X]j=o^i ^ ^' ^^'^ ascer- 
tain that the termwise errors in approximating the Taylor polynomial TQ[d^^) by satisfy 
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analogously as in (37)] 



(39) 



{t - 5.065)^ 



|d(^+4)(5.065) 



■0.065^' <5j (j=0,...,6). 



That the termwise error ( 139|) would not exceed (5j will be guaranteed by A^^- step quadrature 
approximation of the two integrals in (|3]) defining d^^^^\5.065) with prescribed error rjj each. 
Therefore, we set rjj := 6jj\/{2 ■ 0.065-'), and note that in order to have fl5^ 



iv, > n; 



I /T'-f^ II 

l-"5.065,i+4,±lloo 



l-^5.065,j+4,±lloo2 ■ 0.065-'' 



--J'--.?- \; 60-21%- y 60 ■ 2iOj!(5j- 

suffices by the quadrature formula (fTS]) . That is, we must estimate ||-f^5.065,j+4,±lloo for 
j = 0, . . . , 6 and thus find appropriate values of N*. 



Lemma 14. For j = 0, ... ,6 we have the numerical estimates of Table U\ for the values of 
Setting 6j for j = 0, ... ,6 as is given in the table, the quadrature of order 



II H^^ I 

ll-"5.065,j,±l 

500 



Nj > Nj with the listed values of N* yield the approximate values dj as listed in 
TableUl admitting the error estimates (1391) for j = 0,...,6. Furthermore, we have for 
the Lagrange remainder term ||i?6('^^^5 ''^)||oo < 0.0009 =: Si and thus with the approximate 
Taylor polynomial P^it) defined in the approximation \d^^ (t) — P^it)] < 6 := 0.187 holds 
uniformly in [5,5.13]. 



Table 1. Estimates for values of H-f^s.oes j+4,±lloo5 



5„ at; 



and dj for j = 0, . . . , 6. 



J 


\\H^^ II 

ll-"5.065,j+4,±lloo 


6, 




dj 





9.28687 


1014 


0.15 


m 


0.381737508 


1 


2.52880 


10^5 


0.03 


460 


-2.087768122 


2 


6.81644 


10^5 


0.005 


392 


-23.85760346 


3 


1.82039 


10^6 


0.0005 


342 


-140.6261273 


4 


4.82014 


10^6 


0.0002 


196 


-641.9545799 


5 


1.28469 




0.0002 


85 


-2521.387336 


6 


3.80117 




0.0002 


36 


-8940.14559 



Proof. We start with the numerical upper estimation of H^^q^ ^ _^{x) for x E T. For that, 
now we substitute t = 5.065 in the general formula f|T5|) . This results in 

l^5.065,,,±(^)l < 959, 512, 576 ■ v'-''''[j{j - l)(j - 2)(j - 3)f-' + 14.26j(j - l)(j - 2)^-^ 
+ 73.75535j(j - + 163.4085485jF-i + 130.313837600625£^} 

+ 1, 263, 820, 800 ■ t;^-°''^|j(j - - 2)^-''^^ + 12.195j(j - 1)^'"^ + 48.572675jF' 
(41) + 63.105974625^} + 11, 600, 000 ■ t;^"^^ {jf'^ + 5.065F} 

+ 335, 840, 000 ■ t;^-°^^|j(j - 1)F"^ + 8.l3j&-^ + 20.589225^}. 

Otherwise, almost all the functions v''£"^ (with i := |logf|) occurring here satisfy that 
their maximum on < f < 9 is achieved at the right endpoint v = 9. By [HI Lemma 
6], equation (14) this is the case whenever m/s < I/ctq; note that here we consider the 
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degree 6 Taylor polynomial of d^^ , which entails m < 10 while the minimal occurring 
value of s is s = 1.065. So checking the condition m/s < l/cxo ~ 1/0.126 ~ 7.9365.., we 
obtain that f = 9 remains the actual maximum place except for s = 1.065 and m = 9 or 
10. It occurs two times: when max^g[o,9] = (9/(e ■ 1.065))^ = 27126.00128... when 

j = 9, and when j = 10 as power j — 1; together with for s = 1.065 and m = 10, when 
max^g[o,9] t;i-065£io = (io/(e . 1.065))^° = 241857.246.... 

We collect the resulting numerical estimates of in Table [1] and list the corresponding 

values of N* read from formula (HOj) . Moreover, we list in the table the values of dj, too, as 
furnished by the numerical quadrature formula f[T5]) with step size h = 0.001, i.e. N = Nj = 
500 > A^* (j = 0, . . . , 6) steps. □ 

Lemma 15. We have d^^{t) > for all 5 < t < 5.13. 

Proof. We approximate d^^(t) by the polynomial Peit) constructed in fl38l) as the approx- 
imate value of the order 5 Taylor polynomial of d^^ around to •= 5.065. As the error 
is at most 6 = 0.187, it suffices to show that p{t) := P^it) — 5 > in [5,5.13]. Now 
P6(5.13) = 0.188694031... so p(5) = Pei5.13) -6 = 0.188694031... - 0.187 > 0. 

Moreover, p'{t) = P^^it) = . . ^\.A t - 5M5y-^ and p'{5) = -0.806502699... < 0. 

From the explicit formula of p{t) we consecutively compute also p"{5) = —15.96427771... < 0, 
p'"(5) = -103.8163124... < 0,p(^)(5) = -496.9504606... < 0andp(5)(5) = -1940.277873... < 
0. Finally, we arrive at p^^^{t) = d^ = —8940.14559..., which is constant, so p^^^{t) < for 
all t G M. From the found negative values at 5 it follows consecutively that also p^^\t) < 0, 
p^'^\t) < 0, p"'{t) < 0, p"{t) < and p'{t) < on [5, 5.13]. Therefore, p is decreasing, and as 
p(5.13) > 0, p{t) > on the whole interval 5 < t < 5.13. □ 

Next we set forth proving that d'{t) > for all t G [5.13, 5.72]. In this interval we use the 
refined process, applying fl2U]) . Still, for the entire interval we obtain step numbers ~ 550, 
so in order to push down A^ under 500, we divide the interval into 2 parts, and apply 
the method for both sections [5.13,5.33] and [5.33,5.72] separately. That is, we construct 
approximating Taylor polynomials around 5.23 and 5.525. 

So now setting to = 5.23 or to = 5.525, the Taylor approximation of radii tq = 0.1 and 
tq = 0.195, respectively, will have the form 

(42) rf'(t) = V All^t-toy + Rnid',to,t), i?„(rf',to,t):=f— ^(t-to)"+V 

^ j\ (n + 1)! 







Therefore instead of [SI (36)] we can use 

\TD (^1 + +\\ ^ ||gg,«+2,+ ||LM0.1/2] + l|-^g,«+2,-||LM0.1/2] 
-Knl" ,tQ,t}\ ^ — — ■ r, 

(n + 1)! 

/.o\ / |||^g,n+2,+ ||oo + |||gg,n+2 -lloo +i 

^^^^ - (n + 1)! '"^^ 

^ maX|^_tQ|<j.(, ||-ffg,n+2,+ ||oo + niaX|g_tQ|<ro ||i?5,„+2,-||oo n+l 

- (^TfT)! ''^ ■ 

So once again we need to maximize f[T^ . that is functions of the type | logfl^f^, on [0,9]. 
From [SI Lemma 6] and as ^ > 5 for all cases, we obtain for all n + 2 < 5/cro ~ 39.68..., i.e. 
for n < 37 

(44) max \\H.n+2±{x)\\oo< max max logi;r+^ < 9*"+'''' log"+^ 9. 

|5-io|<ro " ^' ' ' |5-to|<roO<^;<9 
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Consider the case to = 5.23. We find (executing numerical tabulation of values for orien- 
tation), that d' is increasing from d'{5.13) ~ 0.0089834... to even more positive values as t 
increases from 5.13 to 5.33. This suggest that it will suffice to approximate d' with an overall 
error just below d'{5.13) ^ 0.0089834.... 

We now chose n = 8, when according to f l44|) ||ifg^io,±(x)||oo < 319,784,241. Therefore 
the Lagrange remainder term fH3|) of the Taylor formula fH2|) with n = 8 can be estimated 
as \R8{d',t)\ < 0.00000176248 ■■■ < 0.000002 =: 69. 

As before, the Taylor coefficients ^'^•''''^^(5.23) cannot be obtained exactly, but only with 
some error, due to the necessity of some kind of numerical integration in the computation of 
the formula (E]). Hence we must set the partial errors 6q, ... ,6s in | '^j-^'-'^^^ '(5-23) (^f_^23y\ < 6j 



such that their sum would satisfy ^j^qSj =: 6 < 0.0089834 in order to have that at least 
d'{5.13) > P8(5.13) — 6 > for the approximate Taylor polynomial 



(45) Ps(t):=^^(^-5.23)^ 



In order to achieve this, we set the partial errors 5q,. . . ,5^ with Y^j=o^j < ^^"^ ascer- 
tain that the termwise errors in approximating the Taylor polynomial Ts{d') by Pg satisfy 
analogously to [8, (37)] 



(46) 



^^^^^»^(t-5.23)^- 



|d(^+i)(5.23) 



■O.V<6, (j = 0. 



We use the refined quadrature (120|1 setting = 500 for all j = 1, ...,8. For this, first we 
need some estimate of the form f|T9l) for |-f^5.2^3j+i,±l all j = 0, . . . , 8. Once such an 

estimate is found with certain exponents {tr,jr) and corresponding coefficients Br,Dr, the 
improved quadrature formula fl2Ul) furnishes an error estimate by means of 



(47) W := W{B, D, t, j) := ^ {P.Q7v(G, t.,>) + Dr Q*^iG, t„>)} • 

r=l 

Namely, the error bound of numerical integration by using our quadrature will then be rjj 
^ and the corresponding termwise error bound becomes 6j — — 



60 -210 -5005' ^ ^ ' 2(j - 1)! ■ O.P-I 

Lemma 16. For j = 0, ... ,8 we have the numerical estimates of Tahle\^ for the values of 
W . Setting 6j as given in the table for j = 0, ... ,8, the approximate quadrature of order 
Nj := N := 500 yield the approximate values dj as listed in Tahle\^ admitting the error 
estimates ( H^ for j = 0,...,8. Furthermore, \\RQ{d' ,t)\\oo < 0.000002 =: and thus with 
the approximate Taylor polynomial Ps{t) defined in fH5]l the approximation \d'{t) — Ps{t)\ < 
6 := 0.004784113 holds uniformly for t e [5.13,5.33]. 
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Table 2. Estimates for values of W and 5,-, dj for j = 0, . . . , 8, with N = 500. 
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J 



1 

2 

3 

4 

5 

6 

7 



estimate 
3.46227 
9.78474 
2.73203 
7.54351 
2.06152 
5.5806 • 
1.4977' 
3.98926 
1.05675 



for W 

•1016 
•1016 
•10^^ 
10^^ 
lO^s 
•10^8 
•10^9 



0.003606534 
0.001019244 
0.000142293 
1.30964 • 10-^ 
8.94756 • 10-^ 
4.84427 • 10-^ 
2.16681 • 10-9 
8.24499 • 10-11 
2.7301 • 10-12 



dj 

0.016265345 
0.084372338 
0.223408446 
-0.41545758 
-8.507038066 
-57.99608037 
-288.5739971 
-1204.823065 
-4474.521416 



Proof. Substituting t = 5.23 in (|T6|) yields 

< Gi-23|G"|5451776{j(j - - 2)(j - 3)| logLl^-^ + 14.92j(j - l)(j - 2)| logL|^-3 

+ 80.9774j(j - 1)1 logLp-2 + 94.465234j| logL^-i + 159.34903641] log L^} 

+ G2-23|G"|7180800|j(j - l)(j - 2)| logL^-^ + 12.69j(j - 1)| logLj^-^ 

+ 52.6787j| logL|^-i + 71.456967] logLpj + G3-23|G"|112000o|j(j - 1)| logLp-^ 

+ 9.46j| logLp-i + 22.1229] logL]^| + G^-23l3872000o|j(j - 1)] logL]^-^ 

+ 9.46j] logL]^-i + 22.1229] logL]^| + G^-2=^11600000 {j\ logL]^-i + 5.23] logL]^} . 

Considering sums of \H^^{C,n)\, this will be estimated by means of Lemma [71 So we insert 

values of j, t and apply Lemma [7] with step number N = 500, getting estimations for W as 

rr. 1 , , , , . W/(60 ■ 210 • 500^) 

IS shown m Table |2j We also calculate Oj = —— — — j — . □ 

2(j — 1)' ■ O.P 

Lemma 17. We have d'{t) > for all 5.13 < t < 5.33. 



Proof. We approximate d'{t) by the polynomial Psit) constructed in fHS]) as the approxi- 
mate value of the order 8 Taylor polynomial of d' around to '■= 5.23. As the error of this 
approximation is at most 6, it suffices to show that p(t) := Psit) — 5 > in [5.13,5.33]. 



Moreover, p'{t) = P^{t) 



;{t - 5.23y-\ Now P8(5.13) = 0.008983405..., and 



(j - 1) 

P8(5.33) = 0.025709673..., so Ps{5.13)-6 > and Ps{5.33)-6 > 0. If we suppose, that p at- 
tain in this interval, that means, the total variation here Var(p) > P8(5.13) +P8(5.33) — 26. 
As Var(p) = J^f.^ \p'\dt, we have an estimation for the integral mean of \p'\ note Ip/ < 
Var(p)/0.2 = 0.12546539.... As it is greater then max(]p'(5.13)], ]p'(5.33)]) and the contin- 
uous function has to attain its integral mean, we have an estimation for total variation of 
p': Var(p') < 2Jp/ — ]p'(5.13) +p'(5.33)]. We also have an estimation for integral mean of 
\p"\: Ipii < Vai(p') /0.2 = 0.43413663.... This process can be continued, till p*-^-* (see Table 
|3]). On the other hand, from the explicit formula of p{t) we consecutively compute also 
p(5)(5.13) < 0, p(6)(5.i3) < 0... Finally, we arrive at p^^^t) = rf8=-4474.521416... How- 
ever, p^^^ is constant, so p^'^\t) < and p^^\t) < in [5.13,5.33]. It means, that p^^\t) 
is decreasing in the interval. It is contradiction, as the calculated lower bound for integral 
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mean of is greater than max(|p*^'^''(5.13)|, |p*^^''(5.13)|), and the function should attain 
this value. □ 



Table 3. Estimates for values of ^'■■'''(5.13), p^^\5.33) and total variation, 
integral mean of p^^^ on interval [5.13, 5.33] for j = 0, . . . , 8 



J 


p(^)(5.13) 


p(-')(5.33) 


Var(p(^)) 


integral mean IpU) 





0.0089834050 


0.025709673 


0.0250930779 




1 


0.061152858 


0.102950595 


0.086827326 


0.12546539 


2 


0.230976823 


0.128352476 


0.508943962 


0.43413663 


3 


0.188714272 


-1.609630427 


3.66852346 


2.544719808 


4 


-3.968140009 


-15.96896377 


16.74813082 


18.3426173 


5 


-34.41704242 


-93.62334897 






6 


-190.4642977 


-431.4289106 






7 


-757.3709229 


-1652.275206 






8 


-4474.521416 


-4474.521416 







In case to = 5.525 numerical tabulation of values gives that d' is positive as t increases 
from 5.33 to 5.72, and rf'(5.33) ^ 0.025709673 . . . , d'(5.72) ^ 0.034577102. We chose n = 9 
then ||ifj_„+2,±(a;)||oo < 1,655,335,712, for this case the Lagrange remainder term (H3i) of 
the Taylor formula (gS]) can be estimated as \Rn{d',t)\ < 0.0000725269 ■ ■ ■ < 0.000073 =: 6g. 
Similarly to ( 145|) and (146|) we now write 

(48) P„(t):=^^(t_ 5.525y, 



7 



(49) 

{t - 5.525)^ 



(5.525) -rfj 



|rf(j+i)(5.525) - dj 



■ 0.195^' <(5j (j = 0, l,...,n). 

We also use the refined quadrature (120|) setting = 500 for all j = 1, ...,9. For this, first 
we need some estimate of the form (IT^ for |-f^5.5^25,i+i,±l ^^^1 for all j = 0, . . . , 9. As before, 
once such an estimate is found with certain exponents {tr,jr) and corresponding coefficients 
Br,Dr, the improved quadrature formula fl20|) furnishes an error estimate by means of W 

W 

defined in flTTI) . with the error of the quadrature being rjj = — -, and the error of 

60 ■ 2^^^ ■ 500^ 

the corresponding term arising from the quadrature becoming 5,- = — — ^ — -. 

^ ^ & 4 5 J 2(j - 1)! ■ 0.195-'-i 

Lemma 18. For j = 0, . . . , 9 we have the numerical estimates of Table [7] for W. Setting 
Sj as given in the table for j = 0, . . . , 9, the approximate quadratures of order N := 500 
yield the approximate values dj as listed in Table^ admitting the error estimates for 
j = 0, . . . ,9. Furthermore, \\Rioid^^ ,t)\\oo < 0.000073 =: 6io and thus with the approximate 
Taylor polynomial Pdit) defined in fl48l) the approximation \d'{t) — P9{t)\ < 6 := 0.0124555 
holds uniformly for t G [5.33,5.72]. 
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Table 4. Estimates for values of W and 6j, dj for j = 0, . . . , 9, with = 500. 



7 
J 


estimate for W 


6, 




dj 





6.89883 


10^5 


0.00715 


6277 


045016622 


1 


1.93082 


10^6 


0.003921976 


0.070827581 


2 


5.34273 


10^6 


0.00105811 


-0.6357179 


3 


1.46288 


10^^ 


0.00015 


8317 


-7.162905157 


4 


3.96656 


10^^ 


2.48926 


• 10"^ 


-45.0748687 


5 


1.06584 


lO^s 


2.60863 


• 10-6 


-220.5767067 


6 


2.84009 


lO^s 


2.2591 ■ 




-922.6394344 


7 


7.50943 


10^8 


1.66398 


■ 10-8 


-3454.236354 


8 


1.97155 


10^9 


1.06486 


■ 10-9 


-11,901.56441 


9 


5.14412 


10^9 


6.01989 


10-11 


-38,448.6079 



Proof. Substituting t = 5.525 in (ITB]) yields 

{H'^^ix)] < G'-'''\G'\5A5l, 776|j(j - - 2)(j - 3)| logL^-^ + 16.1j(j - - 2)| logL^-^ 
+ 94.70375j(j - 1)1 logLp-^ + 120.35253125j| logL|^-i + 222.521187890625] logLp j 
+ G2-525|^/|7^ 180^ 8Qo|^.(^. _ _ 2)1 logLp-^^ + 13.575j(j - 1)1 logL^-^ 
+ 60.426875j| logL^-i + 88.127203125| logLp j + G^-^'^^\G'\1, 120, OOo|j(j - 1)| logLP"^ 
+ 10.05j|logLp-i + 25.000625] logLp j + G3-525^38^ 720, 000|j(j - l)]logL]^-2 
+ 10.05j] logL]^-i + 25.000625] logL]^ } + G^-^^hl, 600, 000 {j\ logL]^-i + 5.525] logL]^} . 

Now the quadrature formula error is to be estimated by means of Lemma [71 So we insert 
the values of j, t and apply Lemma [7] with = 500, obtaining the estimations for W as is 

W^/(60 ■ 2i° ■ 500^) 

shown in Table HI We also get a value of 5j calculating 6j = — — - — . □ 

Lemma 19. We have d'{t) > for all 5.33 < t < 5.72. 

Proof. We approximate d'{t) by the polynomial Pgit) constructed in fH5]) as the approximate 
value of the order 9 Taylor polynomial of d' around to •= 5.525. As the error is at most 6, 
it suffices to show that p{t) := Pgit) — 5 > in [5.33,5.72]. To apply the same method as 
in Lemma fTTl we divide the interval into two parts: [5.33,5,56] and [5.56,5.72]. Moreover, 

p'(t) = P^{t) = ^2%^ . (t - 5.525)-^-i. Now P9(5.33) = 0.025709673..., P9(5.56) = 

U - 1)! 

0.047052108..., and P9(5.72) = 0.034577105... so in these points Pg - 5 > 0. In the interval 
[5.33, 5.56] we get for the integral mean of ]p*^-'^] (estimating with total variation, as in 
Lemma fTT]) . IpU) > max(]p*^''^(5.33)], ]p*^-'^(5.56)]) for j = 1,2 (see Table [5]). The function 
has to attain this estimated integral mean value, so it cannot be monotonic. On the other 
hand p^^\5.33) < for j = 3, . . . , 9. But p^^^ is a constant, that is negative in the entire 
interval, hence p^^^ are also negative in the whole interval for j = 3, . . . , 8. It follows that p" 
is decreasing in the interval, which is a contradiction. 

In case of the interval [5.56, 5.72] the process is similar: Ip/ > max(]p'(5.56)], ]p'(5.72)]), 
and p*^-^^(5.56) < for j = 2, . . . , 9, while p^^^ is a constant, so p^^'^ are also negative in the 
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interval for j = 2, . . . , 9. So p' should be monotonic in the interval, and it is a contradiction. 

□ 



Table 5. Estimates for values of p(^)(5.33), p(^)(5.56), p^^\5.72) and total 
variation, integral mean of p^^^ on [5.33, 5.56] and [5.56, 5.72] for j = 0, . . . , 9 



j 


p(-')(5.33) 


p(-')(5.56) 


p^^\5.72) 


Yar(p^^^) in 
[5.33, 5.56] 


mean /p(j) in 
[5.33, 5.56] 


Yar{p^^^) in 
[5.56, 5.72] 


mean /^(j) in 
[5.56, 5.72] 




1 

2 
3 
4 
5 
6 
7 
8 
9 


0.0257096753 

0.102950466 

0.12835791 

-1.609886701 

-15.96198625 

-93.94395303 

-427.8265683 

-1864.435452 

-4404.08587 

-38,448.6078 


0.047052108 

0.043853873 

-0.915663374 

-8.882443109 

-53.38561447 

-255.0722573 

-1051.102162 

-3894.340881 

-13,247.2656 

-38,448.6078 


0.034577105 

-0.26077396^ 

-3.22675364C 

-21.52543175 

-110.7051556 

-483.1895366 

-1870.009287 

-6506.045571 

-19,399.0429 

-38,448.6078 


0.0478507836 
0.269289432 


0.208046885 
1.170823618 


0.0567182135 


0.354488834 



The last step is to prove that d is concave in the interval [5.72, 6]. 

Lemma 20. We have d"{t) < for 5.72 < t < 6. 

Numerical tabulation of values give that d" is decreasing from (i"(5.72) ^ —0.260774... to 
even more negative values as t increases from 5.72 to 6. In the interval [5.72,6] the second 
derivative of d{t) has the Taylor-approximation 

(50) d"{t) = Y (t - 5.86)^ + RJd", 5.86, t), where 

i=o 

RJd\ 5.86, t) := ^ (t - 5.86)"+' . 

Therefore instead of ^ (36)] we can use 

|i?„(d",5.86,t)| < l|g€.n+3.+ ||LM0,l/2] + l|gg,n+3,-|UM0,l/2] . q^^^.+i 

[n + 1)! 

- (n + 1)! 

^ niaX|g_5.g6|<o.i4 ||-^g,n+3,+ ||oo + maX|g_5.86|<o.i4 ||-^g,n+3,-||oo ^ ]^^n+l 

- (n + 1)! 

So once again we need to maximize ([T^ . that is functions of the type | logf|™f^, on [0,9]. 
From [8, Lemma 6] it follows 

(52) max ||iJg,„+3,±(a:)l|oo < 9^ log"+3 9. 

|5-5.86|<0.14 
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We chose n = 8 then ||iJ^ „_|_3 -|-(x)||oo < 3, 062, 485, 120, for this case the Lagrange remainder 
term §^ of the Taylor iformula dSD]) can be estimated as \Rn{d",t)\ < 0.011209281 ■■ ■ < 
0.00035 =: ^9. 

As before, the Taylor coefficients (ij+2(5.86) cannot be obtained exactly, but only with 
some error, due to the necessity of some kind of numerical integration in the computation of 
the formula Hence we must set the partial errors Sq, . . . , 6s with Y:.=o <^'-= 0.2494, 
say, so that d"{t) < Pnit) + 5 for 

(53) Pn{t) := 5Z ^ (t - 5.86)^' . 



The analogous criteria to [H (37)] now has the form: 



(54) 



rf(^+2)(5.86) -dj- 



it - 5.86)^' 



|rf(-''+2)(5.86) -rfj- 



■•0.14^<5, (j = 0,l,...,n) 



That the termwise error f lM|) would not exceed 5j will be guaranteed by Nj step quadrature 
approximation of the two integrals in ([3]) defining ci^-^+^''(5.86) with prescribed error rfj each. 
Therefore, we set rjj := Sjj\/{2 ■ 0.14-'), and note that in order to have f l54|) 



(55) Nj > n: 



l-^5.86,j+2,±lloo _ 4/ ll-^5.86,i+2,±lloo2 ■ 0.14^'' 



60 ■ 2^% y 60 ■ 2iOj!5j- 



suffices by the integral formula (ITS]) and Lemma 5]. That is, we must estimate ||-f^5j6j+2,±lloo 
for j = 0, . . . , 8 and thus find appropriate values of N*. 



Lemma 21. For j = 0, . . . , 8 we have the numerical estimates of TablelBfor the values of 
ll-^5.86,i.±lloo- Setting 6j as seeing in the table for j = 0, . . . , 8 and 6q = 0.00035, the approx- 
imate quadrature of order Nj > N* with the listed values of N* yield the approximate values 
dj as listed in Table\B, admitting the error estimates fl5^ for j = 0,...,9. Furthermore, 
||-R9((i", t)||oo < 0.01121 =: 6g and thus with the approximate Taylor polynomial Ps{t) defined 
in (153|) the approximation \d"{t) — P8{t)\ < 6 := 0.2494 holds uniformly for t G [5.72,6]. 

Table 6. Estimates for values of H-f^s 86,j+2,±lloo and 6j, N* and for j = 0, . . . , 8. 



J 


II H^"^ II 

ll-"5.86,i+2,±lloo 


5, 


N^j 


dj 





1.07968 


• 10^5 


0.16 


485 


-0.982761617 


1 


2.93801 


10^5 


0.062 


483 


-7.57978318 


2 


7.91604 


W> 


0.015 


453 


-42.74047825 


3 


2.1135 ■ 


10^6 


0.002 


446 


-200.2495965 


4 


5.59555 


10^6 


0.002 


246 


-823.1734963 


5 


1.46998 




0.002 


128 


-3064.925687 


6 


3.83405 




0.002 


64 


-10,561.40925 


7 


9.93361 




0.002 


31 


-34,212.60072 


8 


2.55779 




0.002 


U 


-105,414.5993 
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Proof. We start with the numerical upper estimation of -f^5 86j±(^) ^•'''2 < x < 6. In the 
general formula f|T5|) now we consider the case t = 5.86. 

\H'''{x)\ < 959, 512, 576 ■ v'-''[j{j - l)(j - 2)(j - 3)^"^ + 117.44j(j - l)(j - 2)^-3 

+ 111.5576j(j - 1)^-2 + 309.72742jF"i + 314.40339216^} 

+ 1, 263, 820, 800 ■ - l)(j - 2)&-^ + 14.58j(j - 1)^-''"^ + 69.8588jF~^ 

(56) + 109.931256£^} + 11, 600, 000 ■ v^'^^ {jF"^ + 5.86F} 

+ 335, 840, 000 ■ v^-^^^^j{j - + 10.72jF"^ + 28.4796^}. 

Applying that all the occurring functions of type v'^i"^ have maximum on [0, 9] at the right 
endpoint f = 9 in view of [BJ Lemma 6], we can further estimate substituting i = log 9 
and V = 9. We collect the resulting numerical estimates of ||if^^||oo in Tableland list the 
corresponding values of N* and dj, too, as given by the formulas (153]) and the numerical 
quadrature formula (ITSj) with step size h = 0.001, i.e. N = Nj = 500 steps. □ 

Lemma 22. We have d"{t) < for all 5.72 < t < 6. 

Proof. We approximate d"(t) by the polynomial Ps(t) constructed in (153|) as the approx- 
imate value of the order 8 Taylor polynomial of d" around to '■= 5.86. As the error is 
at most 6, it suffices to show that p{t) := Ps{t) + 5 < in [5.72,6]. _Now P8(5.72) = 

-0.2607741259... so P8(5.72) + 6 <0. Moreover, p'{t) = P^{t) = J^Li , . - 5.86)^-^ 

(j - 1)! 

and p'(5.72) = —3.226759... < 0. From the explicit formula of p{t) we consecutively com- 
pute also /(5.72) = -21.525764... < 0, p"'(5.72) = -110.671188... < 0, p(^)(5.72) = 
-483.626484... < 0, p(^)(5.72) = -1873.40227... < 0, p'-^\5.72)_ = -6804.70822... < 0, 
p(^)(5.72) = -19,454.5568... < 0. Finally, we arrive at p'^^\t) = rf8=-105,414.6... We have 
already checked that ^'^■'•'(5.72) < for j = . . . 7, so in order to conclude p{t) > for 
5.72 < t < 6 it suffices to show p^^^t) < in the given interval. However, p^^^ is constant, 
so p{t) < for all teR. It follows that also p{t) > for all 5.72 < t < 6. □ 
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